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Abstract 

The explanation of the observed galactic magnetic fields may require the exis- 
tence of a primordial magnetic field. Such a field may arise during the early 
cosmological phase transitions, or because of other particle physics related phe- 
nomena in the very early universe reviewed here. The turbulent evolution of the 
initial, randomly fluctuating microscopic field to a large-scale macroscopic field 
can be described in terms of a shell model, which provides an approximation 
to the complete magnetohydrodynamics. The results indicate that there is an 
inverse cascade of magnetic energy whereby the coherence of the magnetic field 
is increased by many orders of magnitude. Cosmological seed fields roughly of 
the order of 10"^*^ G at the scale of protogalaxy, as required by the dynamo 
explanation of galactic magnetic fields, thus seem plausible. 
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1 Introduction 



Apart from the baryon number and the spectrum of energy density fluctuations, the 
physical processes that took place in the very early universe do not have many con- 
sequences that could still be directly detectable today. Most observables have been 
washed away by the thermal bath of the pre-recombination era. One possibility, which 
has recently received increased attention, is offered by the large-scale magnetic fields 
observed in a number of galaxies, in galactic halos, and in clusters of galaxies fl], |^. 
The astrophysical mechanism responsible for the origin of the galactic magnetic fields 
is not understood. Usually one postulates a small seed field, which can then be either 
enhanced by the compression of the protogalaxy, and/or exponentially amplified by 
the turbulent fluid motion as in the dynamo theory . The exciting possibility is that 
the seed field could be truly primordial [^], in which case cosmic magnetic fields could 
provide direct information about the very early universe. 

Early magnetic fields could then play an important role in particle cosmology by 
modifying the dispersion or clustering properties of various particles. One particular 
example is the fate of the neutrino: because of their magnetic moments, Dirac neutri- 
nos propagating in the background of a magnetic field would be subject to a spin flip 
[|], so that a left-handed neutrino can be turned into a right-handed neutrino, giving 
rise to an extra effective neutrino degree of freedom and thereby affecting primordial 
nucleosynthesis. Dark matter particles could also be sensitive to the presence of a 
magnetic field. For instance, axions couple to magnetic fields, but perhaps surpris- 
ingly, it can be shown that despite the coupling, cold axion oscillations are not much 
affected by the presence of a primordial magnetic field 0. 

The issue at hand is then: is it possible that primordial magnetic fields of significant 
strength exist? To answer this, first one has to find a mechanism in the early universe 
which is able to produce a large enough magnetic field. There are various proposals, 
a number of which are based on the early cosmological phase transitions, which are 
discussed in Sect 3. The second problem is to explain how the initial field, which is 
expected to be random as it is created by microphysics and having correlation lengths 
typical to microphysics, can grow up to be coherent enough at large length scales. 
This is a problem in magnetohydrodynamics which is discussed in Sect 4. 

2 Observation and the origin of galactic fields 
2.1 Observing cosmic magnetic fields 

Cosmic magnetic fields can be observed indirectly both at optical and radio wave- 
lengths (for recent reviews, see [Q, 0). Electrons moving in a magnetic field emit 
synchrotron radiation, and both its intensity and polarization can be used to extract 



information about the magnetic field; the radio emission is roughly proportional to 
UeB"^, where Ug is the electron density. However, one first needs to fix the relative 
magnitudes of the electron and magnetic field densities. Usually equipartition of mag- 
netic and plasma energies is assumed, but in some cases free electron densities can be 
determined independently. 

Information about distant magnetic fields, e.g. in clusters of galaxies, has been 
obtained by studying the Faraday rotation of polarized light. The method is based 
on the fact that the plane of polarization of linearly polarized electromagnetic wave 
rotates as it passes through plasma supporting a magnetic field. The rotation angle 
Ax depends on the strength and extension of the magnetic field, the density of plasma, 
and on the wavelength A of radiation. The Faraday rotation measure (RM) out to a 
maximum redshift Zm is defined as 

RM{zm) = ^ = 8.1 X 10^ / ngBi{z){l + z)-^dl{z) rad m~2 , (1) 

where Bi{z) is the strength of the magnetic field along the line of sight in Gauss, ne{z) 
is in cm~^, and 

dl{z) = 10-^Hq\1 + z){l + nz)-^/^dz Mpc (2) 

for a A = Friedmann universe, and Hq is in kms~^Mpc~^ This method requires some 
independent information about the electron density along the line of sight and the 
field reversal scale. In our galaxy this may be deduced from the relative pulse delay 
and frequency of pulsars, but pulsars are too faint to be observed in external galaxies. 
Sometimes ne can be measured independently through X-ray bremsstrahlung, or by 
absorption line transitions in the object to be measured. 

The observed RM is sensitive mainly to a regular magnetic field as small scale 
fiuctuations tend to cancel out in the integral Eq. (|l|). The more tangled the field, 
the smaller RM is. To determine RM accurately one needs to measure at least three 
different wavelengths because the orientation of the polarization plane is ambiguous 
by multiples of vr. The sign of RM allows to distinguish between the two opposite 
directions of the regular field. Magnetic field irregularities may be detected through 
depolarization, which tends to wipe out Faraday rotation at small wavelengths. 

The Zeeman splitting of the spectral lines would provide a direct measure of the 
strength of the magnetic field, but the shifts are very small and this method is appli- 
cable mainly to our own galaxy. It has been used at radio frequencies to study various 
galactic molecular clouds. 

Magnetic fields at the level of few /iG have been detected in galaxies, in galactic 
halos, and in clusters of galaxies. Magnetic field appears to be a normal feature of 
spiral galaxies, where their mean field strengths range from 4/iG in M33 to about 
12/iG 0. Somewhat larger values can be found in the spiral arms. The field in spiral 



galaxies shows often either an axisymmetric or bisymmetric form following the spiral 
structure. It also seems that the magnetic field is not a relative newcomer but existed 
already some 6 x 10^ years ago, as witnessed by the Faraday rotation observed in a 
high (z=0.395) redshift galaxy 0. 

Edge-on galaxies often show vertical dust lines, indicating the existence of vertical 
magnetic field lines escaping from the disc. Synchrotron radiation emitting halos 
confirm this picture. The halo magnetic fields are associated with the outward flow 
of gas and dust ( "galactic wind" ) generated by stellar activity and supernovae in the 

disc g. 

Faraday rotation measurements of radio sources inside and behind galactic clusters, 
combined with electron density estimates obtained from the X-ray spectrum, indicate 
that also many of them have sizeable magnetic fields. An analysis of the RM distribu- 
tion of radio sources seen through or near 50 Abell clusters implies a field strength 
of about 2/xG (Lfrs/10 kpc)~^/^h5Q^, where Lf^s is the typical field reversal scale and h^o 
is the Hubble parameter in units of 50 kms~^Mpc~^; there are indications that Lfrs in 
clusters is actually rather small, of the order of kpc. For example, a recent polarization 
data of a radio galaxy belonging to the Coma cluster implies the existence of a core 

1/2 

magnetic field with B > tangled at scales less than 1 kpc. The fact that 

clusters show a magnetic field which has a strength comparable with the much denser 
interstellar medium of galaxies is somewhat surprising, but would of course be natu- 
ral if the field had a truly primordial origin. Perhaps one piece of evidence pointing 
towards the existence of such a cosmological field is provided by the Hydra A cluster, 
for which the RMs have been determined continuously from within few arcseconds of 
the radio core out to a distance of 40", corresponding to 55 kpc. The high Faraday 



rotations observed |T^, which show strong asymmetry between the north and south 
lobes of the cluster, can be explained in terms of a 6 /iG field coherent over the scale 
of 100 kpc, together with a tangled component of ~ 30 /iG. 

The currently favoured explanation for the origin of the large scale galactic mag- 
netic fields is the a — u dynamo 0, which through turbulence and differential rotation 
amplifies a small frozen-in seed field Bq to the observed fiG field. An initially toroidal 
seed field, which is carried along on the disc of a rotating (spiral) galaxy, is locally 
distorted into a loop by the up- or downward stochastic drift of the gas. As the gas 
moves away from the plane of the disc, the pressure decreases and the gas expands; at 
the same time it is subject to a Coriolis force which will rotate it. The magnetic field 
lines, glued to the gas, will follow and thus a poloidal component perpendicular to Bq 
is generated. The small poloidal loops will reconnect and coalesce to produce a large 
scale field. This is the so-called a-effect. Because the disc does not rotate like a rigid 
body, the field lines will be wrapped and the poloidal field will induce a toroidal field; 
this is the w-effect. Thus the dynamo mechanism does not only produce a large scale 



field but can also, to some extent, predict the shape of the field. 

The dynamo saturates when the growth enters the non-linear regime. A typical 
growth time is of the order of 10^ years, with the rotation period of the galaxy setting 
a lower limit on the growth time. It should however be pointed out that there are some 
indications that the saturation might actually be too fast for a large-scale field to form 



1^. In any cae, the strength of the required initial seed field is rather uncertain, 
but as a rule of thumb one could use a value like 10"^'' G on a comoving scale of the 
protogalaxy (100 kpc). 

Another possibility is that the galactic field results directly from a primordial field, 
which gets compressed when the protogalactic cloud collapses. The primordial field 
strengths needed are however quite large, of the order of 10~^ — 10~^° G. In any 
case, it appears as if a primordial field is required to explain the observed galactic 
magnetic fields, although it is conceivable that a purely astrophysical solution could 
exist as well. For instance, within the nucleus of the "starburst galaxy" M82 there is 
a small but extremely active zone of star formation, which gives rise to considerable 
outflow of material. Associated with it is an outward trasportation of magnetic field of 
about 50 /iG, which manifests itself through field lines connected with the streaming 
motions of the synchrotron emitting relativistic plasma 0, |12[. The field structure 
cannot be attributed to the action of an a — dynamo. This demonstrates that 
other mechanisms exist which are able to amplify and create magnetic fields at least 
in the halos of galaxies. Whether purely astrophysical processes could explain galactic 
magnetic fields remains however to be seen, although some suggestions have been 
put forward. For instance, density and temperature fluctuations could create weak 
magnetic fields via the so-called Bierman battery mechanism which subsequently could 



be enhanced by turbulence |13|. A more exotic suggestion is that rotating black holes 



could be responsible for the seed field [Q. The idea requires asymmetric particle 



emission from the black hole to create a current and hence a magnetic field. 
2.2 Magnetohydrodynamics in curved space 

Before discussing limits on primordial magnetic fields, it is useful to have a brief 
theoretical digression. The behaviour of the early magnetic field is described by the 
magnetohydrodynamic (MHD) equations |T^. To derive them, one assumes that the 
energy-momentum tensor is the sum of an ideal fluid and Maxwell parts which, in the 
absence of viscosity and magnetic diffusivity, reads 

T^" = {P + pWU" + pg^'" + (f^'^F\ - ^(^'^^Fa^F^'^') . (3) 

47r V 4 / 



Here is the four- velocity of the plasma, normalized as U^U^ = —1, and F^jy = 
dfj^Ai, — dyA^ is the electromagnetic field tensor. Note that, as long as diffusion can 



be neglected, the presence of the magnetic field does not change the equation of state 
of an ideal fiuid. The equations of motion for the fiuid arise from energy-momentum 
conservation 

T^",, = -^A.^T^- + T>t^T^^ = . (4) 
The Maxwell equations read 

F^",, = J^, = 0. (5) 

In a fiat, isotropic and homogeneous universe with a Robertson- Walker metric ds^ = 
—dt^ + R^{t)dy?, we may define F^^ in terms of the electric and magnetic fields as 

Fio = RE\ Fij = eijkR^B^, (6) 

where latin letters go from 1 to 3. With this definition the expression for the total 
energy has no -R-factors and takes therefore the familiar form 

T°° = (p + p)7'-p + ^(B2 + E2), (7) 

where 7 = 

The MHD equations take a particularly simple form when using the conformal 
metric ds^ = a?{T){dT'^ — dx^). The Maxwell equations can then be written as 

rfVK 

^ = -VxE, V-B = 0, (8) 

and 

J = VxB-^, V-E = pe (9) 

where the scaled variables are defined as 

pe = a^Pe , B = a^B , J = a^J , and E = a^E. (10) 



Here pe is the charge density. 

The MHD equations in curved space, using conformal metric, thus have a similar 
form to those in fiat space. Let us therefore consider the fiat space for simplicity. The 
evolution of the magnetic field is described by 

rfVK 1 

= V X (v X B) - V X -(V X B) , (11) 
at a 

which follows from the Maxwell equations Eqs. (|],^ and Ohm's law. In a perfect 
conductor a = 00 and field lines are frozen in to the fiuid. If a is finite, the fiuid can 
slip through the field lines. The quantitative measure of the freezing is given by the 



magnetic Reynolds number Rm- It can be defined through the ratio of the two terms 
in Eq. (|llD as 

alV X (v X B)| 
Rm ^ ' ^' ~ avL , (12) 

where L is the length scale and v is the bulk or drift velocity of the fluid. If i?M S> 1, 
the magnetic flux lines will be frozen into the plasma. In the astrophysical context the 
fluid can often be considered as an extremely good conductor with Rm ^ 1. The early 
universe is also a good conductor [|l^, ^ with a ^ T and L ~ Mpi/T"^ so that easily 
Rm ^ 10^°. Large Rm signals the onset of turbulence, which thus could be expected 
to play a role in the early universe magnetohydrodynamics, as will be discussed in 
Sect. 4. 

2.3 Cosmologial fields: limits and observations 

Observing cosmological, intergalactic magnetic fields would naturally be of great im- 
portance as it would strengthen the case for their truly primordial origin. While the 
fields in the neighbourhood of galaxies or clusters of galaxies could conceivably orig- 
inate in the ejected magnetoplasma, magnetic fields in the voids, which are devoid 
of baryonic matter, would be much harder to explain in terms of conventional astro- 
physics. 

Energetic 7-rays can traverse the intergalactic voids, and it has been claimed by 



Plaga |jT8| that the arrival times of extragalactic 7-rays could be used to detect fields 
as weak as 10^^^ G. The idea is that cosmic rays, originating from far-away objects 
such as QSOs or gamma-ray bursters, scatter off the background cosmic magnetic 
field. This gives rise to pair production and a delayed 7-ray which could then be 
observed, and the ratio of prompt to delayed 7's provides a measure of the strength 
of the intergalactic magnetic field. The energy spectrum of 7-rays is likewise affected 



by the intergalactic magnetic field Ultra-high cosmic ray protons would also be 
deflected by the intergalactic magnetic fields so that their arrival times could be used 
pm to set bounds as low as lO^-'^^G. 

A constraint on the strength of a cosmological magnetic field has been obtained by 
considering the rotation measure of a sample of 309 galaxies and quasars, which yields 
the limit |^ B < 10~^G (fiiGhioo/O.Ol)"^, where Qjg is the fraction of the ionized 
gas density of the critical density in the intergalactic medium. The assumption here is 
that the field is coherent over the horizon scale, and the limit is weakened if the field 
is tangled at smaller scales, as in fact would seem hkely. It has also been suggested 
that the power spectrum of the cosmological field could be determined by studying 



the correlations in the RM for extragalactic sources . 

If a primordial magnetic field is present at the time of recombination, it will give 
rise to anisotropic pressures. These would distort the microwave background. Barrow, 



Ferreira and Silk |^ have recently considered the effect of a magnetic field on the 
evolution of shear anisotropy a a in a general anistropic flat universe. For a fluid with 
an equation of state given by p = (7 — l)p (with < 7 < 4/3) it was found that 

H 2--f \ p p J 

where Pga is the gas anistropy, pb = B'^/Sn is the magnetic energy density, and S 
is a constant. The angular anisotropy AT is then directly proportional to aA/H. 
Barrow, Ferreira and Silk then compared the result with the 4-year COBE data set, 
and assuming that the whole observed anisotropy is due to magnetic stresses only, 
concluded that, if the field can be taken homogeneous, B < 3.4 x lO^^G (^ohgQ)^/^. 



This is a more stringent limit than what can be obtained [0, |2^, |26[ from primordial 
nucleosynthesis considerations. 

The anisotropy limit also indicates that the distortions of the Doppler peaks of the 



microwave background [27| due to a homogeneous magnetic field are unobservable, 
although distortions due to an inhomogeneous magnetic field may still be observable. 
The distortions arise because in the presence of a magnetic field the wave patterns 
generating fluctuations in the energy density change. There are magnetosonic waves, 
which in the absence of a magnetic field would correspond to sound waves, together 
with Alfven waves propagating with velocity f a, which however do not induce energy 
density fluctuations but affect the velocities of baryons and photons |^ . The fast 



magnetosonic waves effectively change the sound velocity vs in the baryon fluid by 
t;| — >■ f I + v\ siv? 9, where 9 is the angle between B and the wave vector k, thereby 
changing the Doppler peak structure. 

The Faraday rotation in the polarization of the microwave background could still 
be observable, as has been argued by Kosowsky and Loeb [Q. They estimate that for 
a primordial field corresponding to the present value of 10~^G, coherent on the scale 
of the last scattering surface, there should be a rotation angle of 1° at 30 GHz. Such 
a signal requires a beam size of order 10' and could conceivably be detected in the 
future microwave background maps. 

2.4 Limits from primordial nucleosynthesis 

Primordial nucleosynthesis offers a laboratory for many particle physics ideas. It 
also provides a testing gorund for primordial magnetic fields, which would modify 



the calculated element abundances through two effects p^ , |26| . The extra magnetic 
field energy density would accelerate the expansion rate of the universe, and hence 
increase the predicted abundance of ^He. If one defines the magnetic field at the end 
of nucleosynthesis (T = 10^ K) in terms of the critical field Be = eB/ml = 4.4 x lO^^G 
with B = 7-Bc, one finds numerically that 7 = 3 x 10~^ is sufficient to push 



the relative ^He-abundance above 0.24. There is also a small effect on the ^Li and 
(D +^ He) abundances, but these are negligible. 

However, magnetic fields affect the Hubble rate yet in another, more subtle way 



2J, The statistical distributions of relativistic electrons and positrons change 
because of the background magnetic field. The particles are no longer plane waves but 
occupy Landau levels, with the dispersion relation given by (for B along z-axis) 

E = {pl + ml + 2eBns) ^''^ + m^K , (14) 

where = n + l/2 — Sz with = ±1/2 and n the Landau level, and n is the anomalous 
magnetic moment term, which for fields less than about lO^^G reads k, = eaB /{Airml); 
for subcritical fields, its effect is small. The number density of states in the interval 
dpz is then 

eB 

(2-5n.,o)^rfp. • (15) 

If eB <^ T^, as is natural, a fraction of the electron-positron pairs condense in the 
lowest Landau level. As a consequence, the number and energy densities of electrons 
and positrons increase with respect to the zero field case. 

A primodial magnetic field affects primordial nucleosynthesis also by increasing 



n ^ p reaction rates p4|, |2^ . Hence ^He abundance tends to decrease, partly counter- 
acting the effect due to change in the Hubble expansion rate. It turns out, however, 
that this effect is also small. More important is the modification of the phase space of 
the electrons: increasing number and energy densities causes a decrease in all of the 
weak interaction rates. Taking all these effects into account and requiring that the 



''He abundace should not exceed 0.245, Grasso and Rubinstein obtain the upper 
limit 



B < 1 X lO^^G (16) 



at T = 10 K. A similar conclusion has also been obtained in |]2^. Here the magnetic 
field was assumed to be uniform; if inhomogeneities on the scale much less than the 
horizon are present, as seems likely, the effect of the magnetic field on weak rates could 
be averaged out. In this case Grasso and Rubinstein obtain the limit 

B <lx lO^^G . (17) 
This assumes a magnetic field coherence length Lq in the range 10 ^ Lq ^ 10^^ cm. 



3 Generating primordial magnetic fields 



3.1 Early phase transitions 

There are a number of proposals for the origin of magnetic fields in the early universe. 
Fluctuations in the electromagnetic field in a relativistic plasma are by themselves 
sufficient for generating a small scale random magnetic field To obtain a seed 



field for galactic magnetic fields, one however needs fields with much larger coherence 
lengths. A natural place to look for large scale magnetic fields would be the early 
cosmological phase transitions, such as the electroweak or the QCD phase transition. 

Magnetic fields may arise in a electrically neutral plasma if local charge separation 
happens to take place, thus creating a local current. It has been proposed that this 
could occur during a first order QCD or EW ^ phase transition, which pro- 



ceed by nucleating bubbles of the new phase in the background of the old phase. There 
one finds net baryon number gradients at the phase boundaries, providing the basis for 
charge separation, and the seed fields arise through instabilities in the fluid flow. For 
that one has to require that the growing modes are not damped. Turbulent flow near 
the walls of the bubbles is then expected to amplify and freeze the transient seed field. 
The various hydrodynamical features have been studied in linear perturbation theory 



by Sigl, Olinto and Jedamzik |32], who argued that on a 10 Mpc comoving scale, field 
strengths of the order of 10^^^ G for EW and 10^^'^ G for QCD could be obtained. 
These might further be enhanced by several orders of magnitude by hydromagnetic 



turbulence |3^, |3^, which will be discussed in Section 4.3. 



3.2 Bubble collisions 

In a first order phase transition the phases of the complex order parameter $ = 
pe*®/A/2 of the nucleated bubbles are uncorrelated. When the bubbles collide, there 
arises a phase gradient which acts as a source for gauge fields. The phase itself is not 



a gauge invariant concept, but it has been pointed out by Kibble and Vilenkin |^ 
that a gauge invariant phase difference can be defined in terms of an integral over the 
gradient D^Q. 

Magnetic field generation in the collision of phase transition bubbles has been 
considered in the abelian Higgs model 0. One assumes that inside the bubble 
the radial mode p settles rapidly to its equilibrium value r] and can thus be treated as 
a constant. The dynamical variables are thus G and the gauge field A^. The starting 
point is the U(l)-symmetric lagrangian 

£ = + D^MD^'^y + ^(I'^'l), (18) 



where the potential V is assumed to have minima at p = and p = rj. The simplest 
case is that two spherical bubbles nucleate, one at {x,y,z,t) = (0,0,2o,0) and the 
other at (x, y, z, t) = (0, 0, —zq, Iq), and keep expanding with the velocity v even after 



colliding. Because of the symmetry of the problem, the solutions to the equations 



of motion are functions of z and r = v^t^ — — only, and the U(l) gauge field 
= x^f{T,z). One finds that the generated magnetic field is rapidly oscillating 
and orthogonal to the z-axis and in this case confined inside the intersection region. It 
has a ring-like shape in the (x,y)-plane. It can also be shown that subsequent collisions 
of the bubbles, which now may have a magnetic field inside the bubbles, nevertheless 



lead to a qualitatively similar outcome |36 



There are however two additional important ingredients which need to be taken 
into account: the high but finite conductivity |p!6|, O] of the primordial plasma and the 



fact that in the electroweak phase transition the bubbles will in fact intersect with non- 



relativistic velocities Finite conductivity gives rise to diffusion, the consequence 
of which is to smooth out the rapid oscillations of B, whereas low v will permit the 
magnetic flux to escape the intersection region and penetrate the colliding bubbles, 
where its evolution will be governed by usual magnetohydrodynamics. 

The strength of the generated magnetic field depends on the bubble wall velocity in 
an essential way In the electroweak case the initial growth of the bubble wall is by 
subsonic deflagration, with velocities of the order of 0.05c, depending on the assumed 



friction strength p^. The wall is preceded by a shock front, which may collide with 



the other bubbles. This results in reheating, and oscillations of the bubble radii, but 
eventually a phase equilibrium is attained. The ensuing bubble growth is very slow 
and takes place because of the expansion of the universe. Because the universe has 
been reheated back to Tc, no new bubbles are formed during the slow growth phase. 
Assuming that the abelian Higgs model results are apphcable to the electroweak case, 



one may estimate that [BH] 



5 ~ 2.0 X lO^Y^^ + G, (19) 

where it has been assumed 6o = 1, = e?7 = 100 GeV; the average distance be- 
tween the nucleation centers rave = 9.5 x lO^'^tj^ and the velocity v = 1.2 x 10~^ 
were taken as reference values. Folding in the spectrum of separation of the adjacent 
shocked spherical bubbles averaging over all possible inclinations of the ring-like 



magnetic field, and taking into account the enhancement of magnetic energy due to 
an inverse cascade (see next Section), one arrives at the estimate Brms — 10~^^G for 
the cosmological magnetic field at the scale of 10 Mpc today [p6| . 

There are some indications that the behaviour of the abelian Higgs model does not 



differ qualitatively from the full SU{2) x U{1) case [^H], However, it has also been 



argued by Tornkvist that in the Standard Model the non-abelian nature of the 
gauge fields in fact forbids the formation of magnetic fields in two-bubble collisions, 
which in effect can be gauge rotated away. The result is valid for a specific intial 
bubble configuration and applies only for the initial evolution. The question is thus 



still open, but it seems likely that in bubble collisions magnetic fields are generated 
sooner or later. For example, in a collision of three bubbles magnetic field can no 
longer be rotated away PT . 



3.3 Fluctuating Higgs gradients 

It has been pointed out by Vachaspati that fluctuating Higgs field gradients will 
induce a magnetic field. Such local fluctuations are naturally present at the EW phase 
transition, no matter what its order, since the embedding of the electromagnetic field 
in SU{2) ® U{1)y involves these gradients: 



= -r{V^V, - V^r]) , y^ = ^\f-^ > (20) 

where is the Higgs field. At the electroweak phase transition the field strength Ffp 
can be expected to be constant and 0{my/Y but varies in a random way over larger 
distances. The vector Vi is also random, of course. Its variation is due to the fact that 
the Higgs field makes a random walk on the vacuum manifold of 0. 

To find out the root-mean-square field at large distance scales, one should make 
use of a proper statistical argument. The question then is, what sort of an average 
should one use? The line average, where the gradient vectors are taken to be the 
basic stochastic variables |44|, gives rise to a scaling behavior B ~ 1/\/L. Because 
measuring the strength of B by Faraday rotation of the emitted light involves the 
average along the line of sight, this averaging method could indeed be appropriate in 
these cases. It has however been argued that for the galactic dynamo, only the volume 



average has dynamical significance This has been shown to imply that, in the 

absence of any turbulent enhancement, even under the most favourable circumstances 
the electroweak phase transition can manage only about 10~^°G on the 1 Mpc scale 

m. 



3.4 Vacuum condensates 

Another, a more exotic possibility for generating primordial magnetic fields is based on 
the observation that, due to quantum fluctuations, the Yang-Mills vacuum is unstable 
in a large enough background magnetic field |^ at zero temperature. 

In the early universe the effective energy picks up thermal corrections from the 
fermion, gauge boson, and Higgs boson loops. The detailed form of the thermal 
correction depends on the actual model, but one may take the cue from the SU(2) 
one-loop calculation, where they are obtained by summing the Boltzmann factors 
exp(— /JE'n) for the oscillator modes 

El=p' + 2gB{n + ^) + 2gBS:, + m'iT), (21) 



where 5*3 = ±1/2 (±1) for fermions (vectors bosons). Eq. ( pTD includes the thermally 
induced mass m{T) ~ gT, corresponding to a ring summation of the relevant diagrams. 
Numerically, the effect of the thermal mass turns out to be very important. At high 
temperature, the leading behaviour is given mainly by the bosonic contributions, and 



thus one arrives at the estimate 47 
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cosh(2x) 
smh(x) 



(22) 



where K^{x) = gBf/{AxT'^) + mlx/{gB) and numerically 5V^ ~ 0.02 {gBY which 
serves only to shift the value of B at the vacuum slightly. Thus B ^ is the state of 
the lowest energy even at high temperature. 

A local fluctuation will then trigger the creation of the a new vacuum with non- 
zero non-abelian magnetic field inside a given particle horizon at scales /i ~ T. The 
Maxwell magnetic field is then just a projection of the non-abelian field. The resulting 
magnetic field is exponentially suppressed and can safely be neglected in the case of 
the electroweak gauge group. In the early universe, however, where possibly a grand 



unified symmetry is valid, the suppression may be less severe It is also attenuated 
by the running of the coupling constant. One can then estimate that 

\ \ 2 
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— ~ 3 X 10 



42 



GUT J 



ait) 



G, 



(23) 



where the reference number is for susy SU(5). The field is large enough to serve as 
the seed field. 

Electroweak magnetic condensates have also recently been considered by Cornwall 



4§] , who has suggested that they would give rise to magnetic fields with a net helicity 



via the generation of electroweak Chern-Simons number. Joyce and Shaposnikov |^ 
have argued that a right-handed electron asymmetry, generated at the GUT scale, 
could give rise to a hypercharge magnetic field By via a Chern-Simons term. By 
is coupled to the baryonic current via the anomaly, and would thus induce baryon 
number fluctations which could survive until primordial nucleosynthesis making 
it inhomogenous. These considerations deserve further study. 



3.5 Structure formation and and magnetic fields 

A natural idea is that primordial fields are induced by the mechanism responsible for 
the structure formation, such as inflation. However, the basic problem with inflation 
with regards to magnetic field generation is that the early universe is a good conductor 
so that, ignoring turbulence, the magnetic flux ~ Ba^ tends to be conserved. To avoid 
this, one needs to break the conformal invariance somehow, as was first suggested 



by Turner and Widrow |16|, who considered couplings to the curvature R such as 



RF^ and RA^, as well as photon-axion couplings. Dilaton coupling of the form e*F^ 
has also been considered [0], and interesting field strengths can be obtained at the 
expense of tuning the dilaton coupling strength. If a phase transition takes place 
during the inflationary period, a sufficiently large magnetic field could be created, 
provided however that the phase transition takes place during the final 5 e-foldings 

Cosmic strings have also been suggested to be responsible for magnetic fields 
the field being generated by the vortical motion inside the wakes of strings. There 
are however some open questions related to the stability of the strings, which might 
be remedied if the strings are superconducting [Q. There also exist speculations 
about superstringy origin of the magnetic field. 



4 From microscopic to macroscopic fields 
4.1 Magnetohydrodynamic turbulence 

Even assuming that a primordial magnetic field is created at some very early epoch, a 
number of issues remain to be worked out before one can say anything definite about 
the role primordial fields in the formation of galactic magnetic fields. At the earliest 
times magnetic fields are generated by particle physics processes with length scales 
typical to particle physics. (It has been shown that such fields are stable against 
thermal fluctuations O.) The remaining question is whether it is at all possible for 
the small scale fluctuations to grow to large scales, and what exactly is the scaling 
behaviour of -Brms or the correlator {B{r + x)B{x)). To study these problems one needs 
to consider the detailed evolution of the magnetic field to account for such issues as 
to what happens when uncorrelated field regions come into contact with each other 
during the course of the expansion of the universe. In general, turbulence is an essential 
feature of such phenomena. These questions can only be answered by considering MHD 



in an expanding universe |]T3|, but a big problem is the very large Reynolds number of 
the early universe. Realistic simulations are simply not possible. However, simulations 
can nevertheless be useful in giving an indication of the possible trend. Starting with 
a random 2d magnetic field configuration in a radiation dominated FRW-universe, one 
finds that the magnetic structures coalesce and there is a gradual emergence of larger 
and larger magnetic scales Such a behaviour is encouraging, indicating that the 



field is not really comovingly frozen. On the other hand, because of the different 
topology, 2d MHD is very different from the 3d MHD, and moreover, the Reynolds 
number in the simulation (about 10) is still quite unrealistic. 

Dimopoulos and Davis have also argued out that when two initially uncorre- 
lated domains come into contact, the field at the interface should untangle with the 



plasma bulk velocity v to avoid the creation of domain walls. They propose that the 
correlation length ^ evolves according to 

where H is the Hubble parameter and the velocity v depends dynamically on B and 
should, in principle, be determined from MHD. Qualitatively, this argument also points 
towards a field not frozen in the plasma. 



4.2 Shell models 

Large Rm raises the possibility of magnetohydrodynamic turbulence, or effects not 
visible in local perturbation theory. In ordinary hydrodynamics many properties of 
turbulence, in particular those related to energy transfer and to the spectral properties 
have been studied successfully using a simple cascade model. This is true not only 
qualitatively, but also quantitatively, which is the reason why the cascade model is 
now much used in nonlinear physics . 



The basic idea is that the interactions due to the nonlinear terms in the MHD equa- 
tions are local in wavenumber space, and in fc-space the quadratic nonlinear terms 
become a convolution. Interactions in fc-space involving triangles with similar side 
lengths have the largest contribution. This has led to the shell model which is formu- 
lated in the space of the modulus of the wave numbers. This space is approximated by 
N shells, where each shell consists of wave numbers with 2" < < 2""*"^ (in the appro- 
priate units). The Fourier transform of the velocity over a length scale = 2") 
is given by the complex quantity and -B„ denotes a similar quantity for the _B-field. 
Furthermore, the convolution is approximated by a sum over the nearest and the next 
nearest neighbours, 

2 

Nn{v,B)= Ci,Vn+iBn+y (25) 

Here v and B have lost their vectorial character, which reflects the fact that this 
model is not supposed to be an approximation of the original equations, but should 
be considered as a toy model that has similar conservation properties as the original 
equations. 

Velocity and magnetic fields are thus represented by scalars at the discrete wave 
numbers fc„ = 2" (n = 1, iV), i.e. increases exponentially. Therefore such a model 
can cover a large range of length scales (typically up to ten orders of magnitude). The 
important conserved quantity is -Etotfl"^, where i?tot = / T^^d?x is the total energy. The 
resulting equations of motion read, in conformal time t, 



4 dVn 



Po^ = Nn{vM. (26) 



Figure 1: Spectra of the magnetic energy at different times. The straight dotted- 
dashed hne gives the initial condition [to = 1), the sohd hne gives the final time 
(t = 3 X 10^), and the dotted curves are for intermediate times (in uniform intervals 
of Alog(t-to) = 0.6). 

i = M„(t;,6), (27) 

where po is the background radiation density, and b = Ba^. The functions A^^ and M„ 
each have three terms (combinations of Vj and bj), the form of which is determined by 



the MHD equations, with three free parameters A, B, and C p3 . 



4.3 Inverse cascade 

The numerical study of the cascade model requires of course that the parameters 
A,B,C are fixed so that the model has the same conservation laws as the full-fledged 
MHD. One ignores the detailed evolution of the density by setting p ~ poa'^ so that 

^pov' + ib^) = const. (28) 



The model can then be solved numerically and the results are shown in Fig 



1, where the transfer of magnetic energy to larger and larger length scales is clearly 
seen. This process, the inverse cascade, is due to the nonlinear terms giving rise to 
mode interactions. The initial magnetic energy spectrum was chosen to be given by 
Enit = 0) ~ with the total magnetic energy equal to po- The number of shells 
was = 30 so that length scales differing by ten orders of magnitude were covered. 
It was found that the integral scale, which measures where most magnetic energy is 
concentrated and which is given by 

lo= j{2Ti/k)EM{k)dkl J EM{k)dk, (29) 



Figure 2: Magnetic and plasma kinetic energy spectrum as a function of the wave 
number k in the cascade model for small {u = 10^^) plasma viscosity (left) and large 
{u = 10^) plasma viscosity (right). The initial spectrum has been taken to be flat, 
in accordance with the large-time behaviour of Fig. 1. The highest time t = 10* 
corresponds to the Hubble time 10^^. 

where EM{k) is the magnetic energy spectrum, increases with the Hubble time ap- 
proximately like t^H^^- 

However, around the time of recombination the photon mean free path is very 
large and photon diffusion becomes very efficient in smoothing out virtually all in- 



homogeneities of the photon-baryon plasma This process is often referred to as 
Silk damping, which corresponds to a kinematic viscosity z/ ~ (in natural units). 
Silk damping may thus destroy the magnetic field, as has been noted by Jedamzik, 
Katalinic and Olinto [Q. Therefore one has to follow numerically the evolution of 
the magnetic and kinetic energy spectra in the presence of kinematic viscosity. The 
results are presented in Fig. 2 and the main point can be summarized as follows: 
in the cascade models magnetic energy is transferred to large length scales even in the 
presence of large viscosity. Here the initial magnetic spectrum was chosen to be flat 
in accordance with the large time behaviour suggested by Fig. 1. 



For a sufficiently large viscosity, the inverse cascade stops. One may estimate |^ 
that this typically takes place close to recombination. The results suggest that in the 
real MHD, inverse cascade is operative and is essentially not affected by Silk damping, 
except very late and perhaps for very weak fields. Thus we may conclude that it is 
unlikely that an equipartition exists in the very early universe. A similar conclusion 



can be drawn in a different, continuous model where the inverse cascade can be found 
analytically in an appropriate scaling regime [34, |6l|] . 



5 Conclusions 

Explaining the galactic magnetic fields in terms of microphysical processes that took 
place when the universe was only ten billionth of a second old is a daunting task, which 
is not made easier by the complicated evolution of the magnetic field as it is twisted 
and tangled by the flow of plasma. It is nevertheless encouraging that mechanisms 
for generating primordial magnetic flelds of suitable size exist, and in particular those 
based on the early cosmological phase transitions discussed in Sect. 3 look promising. 
At the same time the fact that there are so many possibilities tends to underline our 
ignorance of the details of the subsequent evolution of the magnetic fleld. The step 
from microphysics to macroscopic flelds is a difficult one because of the very large 
magnetic Reynolds number of the early universe. However, different considerations, 
both analytic approximations, 2d simulations, as well as the full-fledged shell model 
computations which can account for turbulence, seem to point to the existence of an 
inverse cascade of magnetic energy. Moreover, as discussed in Sect. 4.3, the inverse 
cascade is obtained also in the presence of a large plasma viscosity. Therefore the 
primordial origin of the galactig magnetic flelds is quite possible. 

Much theoretical work remains to be done, though. At the same time it is very 
important that progress is made on the observational front. In particular, measuring 
or setting a stringent limit on the intergalactic fleld, which could be possible in the near 
future as indicated in Sect. 2.3, would provide the testing ground for all theoretical 
scenarios. 
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